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XV. On Systems of Linear Indeterminate Equations ■ and Congruences. By Henry J. 
Stephen Smith, M.A., Fellow and Mathematical Lecturer ofBalliol College, Oxford. 
Communicated by J. J. Sylvester, Esq., F.R.S. 

Eeceived January 17, — Eead January 31, 1861. 



The theory of the solution, in positive or negative integral numbers, of systems of linear 
indeterminate equations, requires the consideration of rectangular matrices, the con- 
stituents of which are integral numbers. It will therefore be convenient to explain the 
meaning which we shall attach to certain phrases and symbols relating to such matrices. 
A matrix containing p constituents in every horizontal row, and q in every vertical 



qxp 
A 



, or (when 



column, is a matrix of the type qXp. We shall employ the symbol 

it is not necessary that the type of the matrix should be indicated in its symbol) the 
simpler symbol |A|| to represent the matrix 

■"•I, 19 ■"•!, 25 •"■!, i 



-"■2, 1) Ai, 2? 



A, 



2,p 



A A A 

ri -q, 15 -"-J, 2) 4>P 

If |A|| and |B| be two matrices of the same type, the equation ||A|=|B| indicates that 
the constituents of |A|| are respectively equal to the constituents of |B| ; whereas the 
equation |A|=|B| will merely express that the determinants of ||A|| are equal to the cor- 
responding determinants of ||B|. The determinants of a matrix are, of course, the deter- 
minants of the greatest square matrices contained in it ; similarly, its minor determinants 
of order i are the determinants of the square matrices of the type i X i that are con- 
tained in it. Matrices of the types nx{m-\-n) and »ix(»+i) are said to be of com- 
plementary types ; if ||A|| and ||B|| be two such matrices, we shall employ the equation 

|A|=|B| 

to express that each determinant of ||A| is equal to that determinant of |jB||, by which it 

A 
is multiplied in the development of the determinant of the square matrix 



B 



When 



and 



are dif- 



m and n are both uneven numbers, the signs of the determinants 

ferent : this occasions a certain ambiguity of sign in the interpretation of the equation 
|A|=|B|, which, however, will occasion no inconvenience. If m=n, the matrices J|A| and 
|jB| are at once of the same, and of complementary types; so that, in this case, the 
equation |A|=|B| may stand for either of two very different sets of equations; but this 
mdccclxi. 2 s 
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also is an imperfection of the notation here employed, which it is sufficient to have 
pointed out. If /<; denote any quantity whatever, it is hardly necessary to state that the 

equality 

|A|=/£x|B| 

implies that the determinants of ||A|| are respectively Jc times the corresponding determi- 
nants of ||B|j. 

Let ||P|| be a square matrix of the type nxn>, and |Q|| a matrix of the type nx (n-\-m) 
(where m>0), we shall understand by the matrix compounded of ||Pj| and |Qj|, the matrix 
||X| of the same type as ||Q[|, the constituents of which are defined by the equation 

and we shall write 

||X|H|P||x||Q|; 

in this equation |jQj| is said to be premultiplied by ||Pj|, and |P|j to be post-multiplied by 
||Q|. This definition will suffice for our present purpose ; as the only case of composi- 
tion which we shall have to consider, is that in which the vertical dimensions of the 
matrices to be compounded are all equal, and in which every premultiplying matrix is 
square, so that if an oblong matrix present itself at all in a series of matrices to be com- 
pounded, it will occupy the last place in the series. 

By the greatest divisor of a matrix we are to understand the greatest common divisor 
of the determinants of the matrix. If the matrix be square, its greatest divisor is, con- 
sequently, the determinant of the matrix. A prime matrix is one of which the greatest 
divisor is unity ; i. e. the determinants of which are relatively prime. A prime square 
matrix (i. e. a matrix of which the determinant is unity) we shall call a unit-matrix. 

In any system of linear equations, whether defective or redundant, or neither, we 
shall understand by the matrix of the system the matrix formed by the coefficients of 
the unknown quantities. If to this matrix we add an additional vertical column, com- 
posed of the absolute terms of the equations, the resulting matrix we shall term (for 
brevity) the augmented matrix of the system. 

Lastly, when we have occasion to consider square matrices, the constituents of which, 
excepting those on the principal diameter, are zero, we shall represent them by symbols 
of the form 

lift. £» Ssf-sJi, 

where q„ q 2 , ...q n are the constituents of the principal diameter. 

Art. 2. If every determinant of the augmented matrix of a redundant system of linear 
equations is equal to zero, while the determinants of the unaugmented matrix are not 
all equal to zero, the system admits of one solution, and one only. And in particular if 
the matrix of the system be a prime matrix, the values of the unknown quantities which 
satisfy the system are integral numbers. For these values may be expressed as fractions 
having for their denominators any one of the determinants of the matrix ; and these 
determinants are relatively prime. 
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Let ]|A|| be a given prime matrix of the type nx{n-\-m), |jK|| a given matrix of the 
same type connected with ||A|| by the equation 

|K|=&X|A|, (1.) 

which implies that Tc is the greatest divisor of j|K|| ; then the symbolic equation 

I|K|!=WX||A||, (2.) 

in which \\k\\ denotes a square matrix of the type n x n, will admit of one solution, and 
one only. 

For, to determine Tc Tt „ k r> 2 , . . . h Ti „, the constituents of the rth horizontal row of ||#||, we 
have the redundant system 

-**-»•, i =z -&-l, i «V, lT-^J, t "V, 2~r • • • "T Ai, i «Y, nl 1 .„ . 

is=l, 2, 3,.. .n+m J 

which is involved in the symbolic equation (2.). The matrix of this system is the prime 
matrix ||A| ; and the determinants of its augmented matrix are all equal to zero ; for, by 
virtue of equation (1.), they are equal to the determinants 

-"-r, 15 -"t, 25 • • • • -**r, n+m 
-"-1, 15 -"-1, 25 • • • • -l»-2, n+m 
-**-2, 15 -"-2, 25 • • • • -l»-2, n+m 



-jx 



-"•m, 15 •**-», 25 • • • • -"-« 



in which two horizontal rows are identical. Thus the system (3.), and consequently the 
equation (2.), admits of one solution, and one only. It is evident that the determinant 
of \\k\\ is Jc. The case in which m=0 is not included in this demonstration; its proof, 
however, presents no difficulty, and may be omitted here. 

A particular case of this theorem (that in which w=2) occurs in the ' Disquisitiones 
Arithmeticse ' (see art. 234 of that work). 

Art. 3. If every determinant of the augmented matrix of a redundant system of linear 
congruences be divisible by the modulus, while the greatest divisor of the unaugmented 
matrix is prime to the modulus, the system is resoluble and admits of only one solu- 
tion. For if the modulus be represented by P X Q X R . • ., P, Q, E . . denoting powers of 
unequal primes, one (at least) of the determinants of the unaugmented matrix is prime 
to P, one (at least) is prime to Q, &c. ; whence it may be inferred that the system is 
resoluble for each of the modules P, Q, R.. ., and admits of only one solution for each 
of them; it is therefore resoluble for their product PxQxR .-.5 and admits of only 
one solution for that modulus. 

Let ||K|| denote (as in the preceding article) a given matrix of the type nX(n-{-m), of 

which Tc is the greatest divisor ; and let it be required to find the complete solution of 

the symbolic equation 

||Kf=#||x||A||, (4.) 

in which ||#|| is a square matrix of which the determinant is k, ||A|| a prime matrix of 

2 s 2 
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the same type as j|K||, and in which the constituents of |A|| and \\k\\ are the unknown 
numbers. 

We shall first obtain a particular solution of this equation, and then show how from 
any particular solution the complete solution may be deduced. 

We may suppose that the constituents of any horizontal row of |K|| admit of no 
common divisor but unity ; for if S„ S 2 , ... \ be the greatest common divisors of the 
constituents of the horizontal rows of j|K||, we find 

|jK|HM 2 ,&3, .A||X||K'||, , (5.) 

|K'|| denoting a matrix the constituents of which are derived from those of ||K| by the 
relation 

so that the solution of equation (4.) depends on the solution of a similar equation for 
the matrix ||K'||, in which the constituents of each horizontal row are relatively prime. 



Let then the matrix 



rx(n-\-m) 
K I 



i. e. the matrix 



K K 

X *-l, 15 J - V l, 25 
-"■2, 15 ->*-2, 25 



K,, re+ 
K 



2, n+m 



[l<r<»], 



(r+l)x{n+m) 
K 



admit of a greatest divisor p. 



-**•>•, 15 -t»-r, 25 •••• -*V, n-i 

be a prime matrix, but let the matrix 

Determine a,, a> 2 , .... u r by the system of congruences, 

K M »,+K 2)i <y 2 +... K M <y r =K,. +1; ,., mod, ^,1 

t=l, 2, 3, ...n+m J ( -) 

(which, as we have just seen, is always resoluble), and in ||K|| replace the constituents 
K r+M by the numbers 

~[K r + 1, i — 2 S= 1 &» s K Sj J ; 

we thus deduce from ||K|| another matrix ||K"|| connected with it by the relation 
|K|=/&x|K"|, and such that the matrix of its first r+1 horizontal rows is prime. By 
proceeding in this manner, we shall at last obtain a prime matrix ||A |j, which satisfies 
the equation |K|=/£X|A |; we may then, by the method of the last article, determine a 
square matrix \\k \\ satisfying the equation 

||K||=|N|X||A |J, (8.) 

and thus obtain a particular solution of the proposed equation (4.). 

To deduce the general solution of that equation, let ||#j|| and j|A,|| be any two matrices 
satisfying it. We have therefore the equality 

||*,||X|[A I |H|£ ||xl|AJ|, (9.) 
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which evidently implies that 

NHA.,1; (10.) 

whence, by the theorem of the last article, 

lA.HMIxINi, (ii.) 

||«|| denoting a unit-matrix. Combining (9.) and (11.), we find 

|^IXHX||A ||=WX|!A ||, ....... (12.) 

whence, by the same theorem, it follows that 

IMx|Ml=#.|; (13.) 

or, which is the same thing, 

INNNIxH" 1 , (H.) 

j]a|j -1 denoting the matrix reciprocal to ||«|. The complete solution of equation (4.) is 
therefore contained in the formula? 

« A INI« llx« A «« 1 . (15 . } 

||K||=||K ||x|| a I- 1 ] 

j]a|| denoting an arbitrary unit-matrix of the type n X n, and |A |, ||# || being any two 
matrices that satisfy the equation. 

In this, as in the preceding article, we have for simplicity excluded the case in which 
m=0, and the matrices |K|| and ||A| are squares. But it is readily seen that no exception 
is presented by this particular case. 

Art. 4. Let 

J- ..... . (lo.) 

«=1, 2, 3, ...n j V ! 

represent a system of indeterminate equations of which the matrix is |A||. We shall 
suppose that the determinants of ||A|| are not all equal to zero, i. e. that the system is 
independent ; so that its index of indeterminateness (or the excess of the number of 
indeterminates above the number of really independent equations) is m. If we take r 
solutions of the system, for example the solutions 

/v* /v* m /i* ~\ 

8=1, 2, 3, ...r | ( ' j 

it is evident that if r >m, the determinants of the matrix |jar|| are all equal to zero. If 
r< m, and if the determinants of the matrix ||#| be not all equal to zero, the solutions 
(17.) are said to form a set of r independent solutions; if r=m, they form a complete 
set of independent solutions. A set of relatively prime solutions is an independent set of 
which the matrix is prime ; a complete set of relatively prime solutions may be called, 
for a reason which will presently appear, a fundamental set of solutions. It is always 
possible, in an infinite number of ways, to assign complete sets of independent solutions 
of a system of equations of the form (16.). Among the methods by which this may be 
accomplished, we shall select one which depends on the following principle : — 
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If 



l fX(m+»') : 
a 



represent any matrix of the type r X (m+r), the determinants of which 
are not all equal to zero, and if X„ X 25 ... K m+r be integers which satisfy the equations 

......... (18.) 



2j.=i a itk y<- k — 0;1 



i-1, 2, 3 ...rj * 
while «,. +1) i, «,.+!, a , «,. + i, 3 .... flV+n.m+r are integers satisfying the inequality 

^*= "V+l, A A* U, ........ . 



(19.) 



the determinants of the matrix 

(r+l)x(m+r) 
a 
are not all equal to zero. 

For if 2*"" ,+r a r+i, k ^—^ it is evident that by combining this equation with the equa- 

rX(m-j-r) 

j(r+l)x(m+^) 

I <z 

minants do not all vanish, neither can the latter. 

Let, then, a m>n a m>2 , ... a m>n+m represent any particular solution (other, of course, 
than that in which every indeterminate is equal to zero) of the system (16.) ; and let 



tions (18.), we may express each of the determinants #X 
linear function of the determinants of 



m succession as a 



If, therefore, the former deter- 



the system 



A, 



+ 1, n+m 



be integral numbers satisfying the inequality 



>, 






< 



A', 1 *1 I ■"% 2 ^2 I • • • "T-^-i, n+m ^n+m — " 

»s=l, 2, 3, ...»+l 



(20.) 
(21.) 



(which is obtained by the addition of a single equation to the system (16.)) is itself an 
independent system, as appears from the principle just enunciated ; its index of inde- 



terminateness is therefore m— 1. Let 



(m— l)x{n-\-m) 



represent a complete set of 



independent solutions of (21.); it may then be inferred, from a second application of 

the same principle, that ^ ' represents a complete set of independent solutions 

of the proposed system (16.). Thus the determination of a complete set of independent 
solutions of a system of which the index of indeterminateness is m, depends on the 
determination of a similar set of solutions for a system of which the index is lower by 
a unit. By successive reductions, therefore, we shall at last arrive at a system of which 
the index of indeterminateness is unity, the complete solution of which is of course 
immediately found by evaluating the determinants of its matrix. 

The practical application of this method supposes only that we can always assign a 
particular solution of a system of the form (16.) or (21.). And this, it may be observed, 
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can always be done, either by trial, or by other obvious and not unsymmetrical expe- 
dients. 

Art. 5. If ||«| represent the matrix of a complete set of independent solutions of the 
proposed system (16.), and \\b\\ be any matrix of the same type as ||a||, and connected with 
||#|| by the equation 

IN=Pfx||4 (22.) 

in which \\k\\ denotes a square matrix of which the determinant is not zero, it is evident 
that the constituents of \\b\\ are also a complete set of independent solutions. And, con- 
versely, if \\b\\ be the matrix of a complete set of independent solutions, ||a|| is also the 
matrix of a similar set. For if |K| be the matrix composed of the first minors of |[#||, so 
that 

||K|xpHM»*.----ll» 

we have from (22.), 

pC||x||5||=#, Jc, ... k, m-x\\a\\; 

from which it appears that \\k, k, . . . ||x|H|, and therefore ||«j| itself, is the matrix of an 
independent set of solutions. 

This observation enables us to obtain a complete set of relatively prime solutions, as 
soon as we have obtained an independent set. If ||5|| be the matrix of the independent 
set, we have only to determine, by the method of art. 3, a square matrix ||#j|, and an 
oblong prime matrix |«|, satisfying the equation 

IH|=#l|x|H|; 

the constituents of |ja| are then the terms of a set of fundamental solutions. 
Or again, if in art. 4 we employ, instead of the inequality (19.), the equation 

2i=i «,-+i, j.Xj=1, (23.) 



it is easily shown that if 



rx(n-\-m) 
a 



(r-\-l)x{n-\-m) 

a 



is also 



be a prime matrix, 

a prime matrix; so that, by following the method of that article, we may obtain 
directly a set of fundamental solutions of any proposed system. Only, it will be 
observed, that in this mode of obtaining such a set, we suppose that we can assign par- 
ticular solutions, not only of systems of the form (16.), but also of equations of the 
form (23.). 

Art. 6. The importance of fundamental sets of solutions in the theory of linear inde- 
terminate equations is evident from the following proposition : — 

" If ||«|| represent a set of fundamental solutions of the system (16.), the complete 
solution of that system is contained in the formula 

1=1, 2, 3, . . . n+m) 

in which £„ £ 2 , . . . | m are absolutely indeterminate integral numbers." 

For it is evident that every set of numbers included in (24.) satisfies (16.) ; and, con- 
versely, if 0Vh,i, a m+h2 , a m+hn+m be any solution of (16.), the determinants of the 
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matrix 



are all zero, while the matrix 



mx(n-\-m) 

a 



is prime ; whence, 



(m-\-l)x(n-\-m) 
a 

by a principle employed in art. 2, the system 

i=l, 2, 3, n-\-m 

is satisfied by one, and only one system of integral values of |„ | 2 , . . . | m ; or, which 
is the same thing, the numbers a m+l>1 , a m+hVS , ... a m+hn+m are included in the for- 
mula (24.). 

It may be added, that no fractional values of |„ | 2 , . . . |,„ can give integral values 
to #„ s 2 , . . . x H+m ; and that the same values of #,, oc 2 , x 3 , ... x u+m cannot arise from 
different values of |„ | 2 , ... | m . 

The converse of the proposition just established is also true ; *'. e. 

" If the formula 

*,=:Cr&«M .......... (24.) 

represent every solution of an indeterminate system of equations, the matrix \\a\\ is a 
prime matrix." 

For if \\b\\ represent a set of fundamental solutions of the indeterminate system, we 
may express the constituents of \\b\\ as linear functions of the constituents of ||a|j, by 
means of the equations (24.), so as to obtain an equation of the form 

W=|5ilx|!4 

|j|jj denoting a square matrix ; whence it immediately appears that \\a\\ is a prime matrix, 
and |'||| a unit-matrix. 

Thus, if we apply Euler's method for the resolution of indeterminate equations to 
the system (16.), we obtain, as the final result of the process, a system of equations of 
the form (24.); and as it is demonstrable, from the nature of the method itself, that 
these final equations contain the complete solution of the proposed system, their matrix 
is a prime matrix. 

If ||a|| and||#|| be any two sets of fundamental solutions cf the same system, we shall 
have the equation 

Pli=flxj|4 

|||] denoting a unit-matrix. The matrices, therefore, of all sets of fundamental solutions 
are deducible, by premultiplication with unit-matrices, from the matrix of any given set 
of such solutions. 

Art. 7. If ||«| and \\b\\ represent two complete sets of independent solutions of the 
same system, the determinants of ||«| and \\b\\ are evidently connected by the relation 
fix\a\=+aX\b\, a and (3 denoting the greatest divisors of \\a\\ and \\b\\ respectively. A 
similar relation subsists between the matrix of the system and the matrix of any com- 
plete set of independent solutions of it. 

Let |A|| and ||«|| represent those matrices respectively, K and k their greatest divisors ; 
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the relation in question is expressed by the formula 

kx\A\=Kx\a\, (25.) 

where it is to be remembered that the types of the matrices ||A|| and ||«|| are comple- 
mentary; so that, as has been already observed (see art. 1), there is an ambiguity of 
sign in the equation (25.). 

To obtain its demonstration, let Q and q denote the sums of the squares of the deter- 

A 

minants of ||A| and ||«| respectively, and consider the determinant 

is certainly not zero, for multiplying it by itself, we find 

A 



a 



This determinant 



Let, then, 
A 



a 



= QXq (26.) 

be multiplied by any determinant of ||A||; for example, by 2+A,,,, A 2>2 , 
Observing that 2+A,, „ A 2>2 , . . . A Bj „ may assume the form 



A;, i , A 2> i 



*->», i 



,0,0, 



Ai, 2 , A 2> 2 , . . . A_ 2 ? 0, 0, 



■&-\,n > As, » , . . . A. n>n , U, U, U 

A-i, n + i , A 2j n+1 , ... A Mj M+] , 1, U, U 

-"■i, n +2 ? A. 2| M+2 , . . . -A. Mj n+2 , U, 1, U 

a a Ann i 



we obtain the equation 
A 



X Z'iZ ■&■!, u A 22 , ... A- n<% — UXiTfli.n+i) ^2, m+25 



a„ 



;, «+»*) 



(27.) 



in which we may permute the second set of indices in any manner consistent with the 



condition that 



should not change its sign ; so that we may write 
A 



(28.) 
A 



X|A|=Qx|a|, 

the correspondence of the determinants in |A| and \a\ being fixed by the matrix 

The equation (25.) is an immediate consequence of this result; and if in that equation 

we suppose the correspondence of the determinants to be still fixed by the matrix J * 
we shall have to write 



a 



or 



according as 
MDCCCLXI. 



&x|A|= KxK 

#X|A|=— Kx|«|, 
is a positive or negative number. 

2 T 
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Art. 8. From the preceding principles we may deduce the solution of the following 
problem, which admits of important applications in other parts of arithmetic : — 

" To find all the matrices of a given type, of which the determinants have given 
values, not all equal to zero." 

Two particular cases of this problem (those in which the matrix is of the type 2x3 
and 2x4) occur in the ' Disquisitiones Arithmeticae' (see arts. 279 and 236). In both 
places Gauss has suppressed the analysis of the problem, and has only given a synthe- 
tical demonstration that its conditions are satisfied by the solution he assigns. This, 
indeed, in art. 279, he expressly observes. He has also suppressed his method of 
deducing the complete solution from any particular solution, — an omission, however, 
which may probably be supplied by a comparison of art. 234 with art. 213, i. The very 
general and most important case, of a matrix of the type nX (n-\-l), has been subse- 
quently treated of by M. Heemite *. 

Let \\x\\ denote a matrix of the type n X (%+m), of which the constituents are abso- 

lutely indeterminate quantities ; writing X for "tt-jj ? we shall represent its determi- 
nants by X,, X 2 , X. K . . If m>l, these determinants are not all independent, but 

are connected by certain identities of the form 

<J>(X I5 X 2 , XJ=0, (29.) 

<E> denoting a rational and integral homogeneous function with numerical coefficients. 
If, therefore, C„ C a , . . . C A be a given set of integral numbers, which can be represented 
as the determinants of a matrix of the type nx(n+m), these numbers will satisfy every 
relation of the form (29.) ; so that the identity 

<D(X„X 3 , ...X,)=0 

will involve also the numerical equation 

<J>(C 1? C 2 , ...C x )=0 (30.) 

To obtain a convenient notation for C„ C 2 , . . . C x , let us imagine that we have formed 

a square matrix of the type (n-\-m) X (n-\-m) by the addition of to horizontal rows to the 

matrix |#|| ; if, in the development of the determinant of this matrix, the coefficient of 

X,- be the determinant 

r=l, 2, 3, ... to 

s=l, 2, 3, ... to 



X n +r,n s 



(fa, fa, . . . (jb m denoting to of the numbers 1, 2, . . . n-\-m), we may represent X,- and C, 
by the symbols \jn, x , fa, . . . fa^\ and (/*,, //» 2 , . . . p m ) respectively ; observing, however, that 
if two of the numbers fa, fa, . . . are equal, the value zero is to be attributed to each of 
these symbols. 

If r denote one of the numbers 1, 2, 3, . . . n, the determinants of the matrix obtained 

* Oeeile's Journal, vol. si. p. 264 ; see also Exsenstein, ibid. vol. xxviii. p. 327. 
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by adding the horizontal row 



to the matrix 
equations of the form 



&r, 15 <Br, 2J • • • • tt r , »+ 

nx(n-\-m) 

x 



, are identically equal to zero, We thus obtain XX 



n+l 



2, = ™ + * p, jk.,, p„ ... i » m _ 1 ] x rt { =0, (31.) 

/;*!, (a 2 , . . . \x, m _ x representing any combination of m— 1 of the numbers 1, 2, 3, . . . m-\-n. 
In connexion with these equations, consider also the similarly formed system, 

2 , i ir + "(i^i ) ^"-^ m -i)^=0 . (32.) 

This system, which is in appearance redundant (containing XX — tt equations, and 

only m-\-n indeterminates), is in reality defective, and is equivalent to m independent 
equations. For if (k x , Jc 2 , ... Jc m ) be one of the given numbers C which is not equal to 
zero, the partial system of m equations 

2 i=1 (l, fC x , fc 2 , ... rCj-i, Kj+u - •• " ; m)yi == 



i=l,2,3,...m j (33 - } 

is certainly an independent system, because the determinant of the coefficients of 
y h , y h , ...y hm is (#„ k 2 ,... k m ) m , and is therefore different from zero. Again, every equa- 
tion of (32.) which is not already comprised in (33.), may be obtained by linearly com- 
bining the equations of that partial system. To verify this assertion, let 

2[», (*» p*, ■ ■ ■ /*„-,>,, i=0 "1 - ■ 

r=l,2,8,...» J (84) 

be the system of n equations obtained by attributing to r the n values of which it is 
susceptible in any one of the equations (31.). Eliminating from this system those n— 1 
determinants [i, /»„ f* s , ... p m _i] in which i has a value not included in a set of m+1 
numbers v,, c 2 , . . . c m+1 , arbitrarily selected from the series 1, 2, 3, . . . n-\-m, we obtain a 
relation, which may be expressed in the form 

2;.:r +1 (-iy[wi,p 2 , •••f* OT -.]x[> 1 , »„ ...»,.„»,+„ ...wi]=o, • • (35.) 

mn>? . . . 

representing ( m+ n( w + i\ equations, since the sets 

f*n f*a? • • • pm-ii 
'n v 2 ■> •' • v m+\ 

may respectively denote any sets of m— 1 and m+1 numbers taken from the series 
1, 2, . . . m-\-n. Since (35.) is of the form ^(Xj, X 2 , . . . X x )=0, we may at once infer 
the corresponding relation, 

2;ir +1 ( — iy(v { , (*„[**, .../*„,_,) X (»»,»» ••• »i-i,»i+i,.»» m+ i)=0, . . (36.) 

by means of which any one of the equations (32.) may be deduced from the equations of 
the partial system (33.). Thus, if we multiply the equations of that system taken in 

2 t2 
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order, by the determinants (—iy(kj, h u h 2 , . . . h m _ x ), and add the results, we obtain 






(i,h l ,h n ,..h m _i)y i =0, 



i. e 



(k s , k 2 , ... k m )~2 { , 
. since (k„ fc 2 ,... fc m ) is not zero, 



Let 



represent the matrix of (33.), or of any other independent system 



The system (32.) is therefore equivalent to a system of m independent equations. 
m x (m-\-n) 
7 
equivalent to (32.) (the determinants of all such matrices are proportional); let 

r\, r 2 , . . . I\ be the determinants of ||y| ; |||| and H l5 E 2 , . . . E x the matrix and deter- 
minants of the system similarly derived from (31.). By the theorem of art. 7, we have 



x 



x|£=2E 2 xH, 



(37.) 



or observing that 



x 



:S.EX, and that (37) is an identity of the form 4>=0, 

2.rcx|r|=2Px|C|, (38.) 

where |C| symbolizes the numbers C„ C 2 , .... C A , which correspond to the determinants of 
|y|| in the same inverse order in which in equation (37.) the determinants of \x\ corre- 



spond to those of ||||. But if 
of (32.) or (33.), we have also 



wX(wi-f-%) 

9 



represent a system of fundamental solutions 



X|y|=2Px|0., • • f (39.) 



whence, combining (38.) and (39.), and representing the greatest common divisor of 
C„ C 2 , ... C A by c, we find 

cx|fl|=|q (40.) 

If, then, ||c|| denote any square matrix of determinant c, and of the type n X n, the formula 
\\c\\ X P\\ contains the complete solution of the problem. 

If y represent the greatest divisor of |y|, we infer from (38.) 

<?x|y|=yX|C|, (41.) 

whence, if ||y'|J be a prime matrix of the type mx(m-{-n) satisfying the equation 

|y|=yx|y'| (see art. 3), 
we find 

|C|=cx|y'| (42.) 

The preceding analysis enables us therefore to obtain simultaneously the representation 
of the determinants |C| as the determinants of two complementary matrices, of the types 
n X (m+») and m X (m-\-n) respectively. We have thus two distinct methods of arriving 
at the solution of the problem, of which one requires the determination of a set of 
fundamental solutions of a system of linear equations ; the other the reduction (by the 
method of art. 3) of a given matrix to a prime matrix. The greatest divisor of ||yj|, 
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which we have represented by y, is evidently (k u k 2 ,... Jc m ) m ~ l X c. If, therefore, C, one 
of the given numbers C n C 2 , ... Cy, be a unit, we have only to take C for (k t , k 2 ... k m ), 
and we shall immediately obtain a matrix |y| of the type mx(w-f-w) satisfying the 
equation 

And similarly might a matrix of the type n X (m+n), satisfying the same equation, be 
written down without calculation. 

Art. 9. The importance of the case, in which m=l, is so great, that we maybe 
allowed to point out the identity of the solution obtained by the preceding method 
with that already given by M. Hermite. Let, then, C„ C 2 , ... C„ +l represent the deter- 
minants of a matrix of the type »x(»+l) taken in their natural order («*. e. so taken 
that if the matrix be completed by an additional row of constituents, 

Cj, C^ ... C M +j, 

the value of its determinant would be 

CA+CA+C3C3+ ...c n+1 C n+1 ). 

We have then to obtain a set of fundamental solutions of the equation 

0^+0^+0^+... +C M+1 ^ +l =Q (43.) 

Such a set may always be obtained by the following particular method. Supposing 
that 0, is not zero, consider the equations 

0=0,^+0^ 



1+\ 



0=O 1 y 1 +C 2 y 2 +C 3 ^+ +G n+l y n , 

and take a particular solution of each of them, assigning to the last indeterminate in 
each, the least value (zero excepted) of which it is susceptible. If we denote by A* the 
greatest common divisor of C„ C a , . . . C A , so that A,=C„ A n+1 =c, it is evident that the 
value of y k+l in the equation 

Cj^+0 3 y 2 + . . . +C A+1 y ft+1 =0 

and if in the same equation we represent the values of y lt y^ ... y k by 



will be 



A* 



the matrix 



**, 15 ^k, 25 **, 35 ' ' > ' 1"k, *? 



r ~ 



r 2»»' r 2.2» A3 ' 



' 3, 15 '3,2? '3,35 ^ 



**«, U *V 25 5*», 35 *», 4 











Vh 



(45.) 



306 ME. H. J. S. SMITH ON SYSTEMS OF LINEAR 
will represent a set of fundamental solutions of (43.). For, in the first place, it repre- 
sents a set of independent solutions ; because its first determinant is — ' X -r 2 X — > 

1 A 2 A 3 A„ +1 

A C 

or t— K or ~~ ! ' therefore its determinants are proportional to C„ C 2 , . . . &c. ; or since 

the first of them is — , they are respectively equal to the numbers 

Cj C 2 L'g t^m + I 

' 9 9 9 * . • 9 

c c c c 

which admit of no common divisor. 

To obtain a set of values for the constituents r itj , which occur in the matrix (45.), we 
may form the series of equations 

\ C 3 +^ 2 A 2 =A 3 



*- ra -iC M +p„A n = A a+ , ^ 



(46.) 



It will then be found that the equations (44.) are satisfied by the values of r comprised 
in the formula 

ri,j=-h-i to- •••,%-! aE! U^f]'> (47.) 

and on substituting these values in the matrix (45.), it will coincide, after an unim- 
portant modification, with that occurring in M. Hermite's solution of the problem. 

But, in practice, the simplest method of obtaining a solution of the problem con- 
sidered in this article, is to solve the equation (43.) by Euler's method, and to employ 
in the place of the matrix (45.), the matrix of the set of fundamental solutions thus 
obtained (see art 6).j 

Art. 10. Another problem, closely connected with the preceding, and of no less fre- 
quent application, has also been completely solved by M. Hermite*; but as it may serve 
to illustrate the utility of the methods employed in this paper, we shall venture to resume 
and generalize it here. The problem is 

" Given a set of n-\-\ numbers C„ C 2 , . . . C n+l without any common divisor, to assign 
all the matrices ||#|| of the type nx(n-\-l) which satisfy the equation 

x ~~ 

Let c„ c 2 , . . . c„ +1 be any particular solution of the equation 

C 1 y 1 +C 2 y 2 +....+C M+] y B+1 =l (48.) 

(which is always possible because C„ C 2 , . . C„ +1 are relatively prime); and let §yj| repre- 
sent a set of fundamental solutions of the equation 

c«yi+<^«+.-.+<>.+iy, + i=Q . (49.) 

* Liouville, vol. xiv. p. 21. 
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Then, if ||w| represent any unitanatrix of the type nxn, and X„ X 2 , . . . K n absolutely inde- 
terminate integers, the complete solution of the problem is contained in the formula 

IM|x||n/+M^| ' 

*=1, 2, 3 ... » 1 , (50.) 

j=l,2,S...n+l 
For if ||a?|| be any one of the matrices contained in that formula, it is readily seen that 



=1, M is included in the 



=C 1 o,+(Vi+ . . +C„ +1 <Vh=1. 

C 
Conversely, if Ml be a matrix satisfying the equation : 

formula (50.). To show this, we observe that the complete solution of equation (48.) is 
contained in the formula 

y j =c J +?= n J u \ i ,j=l,2,...n+l, (51.) 

in which ||4|| is any set of fundamental solutions of the equation 

C^+C 2 y 2 +....+C a+ #„ +1 ==0, (52.) 

and X„ X 2 , . . X„ are indeterminate integers. The complete solution of the same equa- 
tion (48.) is therefore supplied by the determinants of the matrix ly^.-f X.-CJ. For those 
determinants may be represented by the formula 

<H-2;:?ftjT!W==l, 2, 3, . . .n+l, 

CI 



in which [i,f\ symbolizes a first minor of the determinant 

,C 



, so that 



d 



ujy 



'tyt.i 



But the numbers [i, 1], [i, 2], . . . [^ n-\-l~] satisfy (52.) for every value of i; and, since 
C 



:1, the determinants of the matrix 



, |r . t=l,2, 3, ..» 

L^l- TOO „J_1 

j=l, 2, 3, ... w+1 



(53.) 



a matrix satisfying the equation 
values can be assigned to X„ X a , 



are the numbers C„ C 2 , . . . C n+1 , and are therefore relatively prime. It follows from 
this that (53.) represents a set of fundamental solutions of (52.) ; i. e. that the com- 
plete solution of (48.) is represented by the determinants of fly,-, ; -+*iCJ- If then j|#|| be 

C 

=1, since the determinants of ||#|| evidently satisfy (48.), 

.. . X„ which shall verify the equation 
whence it follows that 

WI=IMIxIny+ x Al' 

\u\ denoting a unit-matrix, i. e. |#|| is one of the matrices included in the formula (50.). 
The result incidentally obtained in the foregoing analysis, that the complete solution 
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of an equation of the form 

L> 1 X 1 -\-L'^t! !l -\- -\-L> n+1 W n+1 =:i 

can be exhibited in the determinantal form (50.), is occasionally useful. 

The preceding problem is a particular case of the following more general enunciation : — 
" Given a prime matrix ||C|j of the type mx(in-\-n), to find all the matrices \\x\\ of the 

type nx(in-i-n) which satisfy the equation 

C 



x 



= 1." 



(54.) 



Let |y|| be a matrix which satisfies (54.), let the numbers frj represent absolute inde- 
terminates, and \\u\\ any unit-matrix ; the complete solution of the problem is contained 
in the formula 

N=IMIxfir+ 2 f* c ll> ......... 



(55.; 



where ||y+2^C|| represents the matrix, 



i=l, 2, 3, ... n 
j=l, 2, 3, ... n+m. 

For if \\x\\ be a matrix satisfying the equation (54.), we have 



c 

X 


7 




c 

X 


=Hx 


c 

7 



and consequently 



||t| denoting a unit of the type (m-\-n) X (m+w). But because the first m horizontal rows 



m 



and 



are identical, it is evident that 



i=l, 2, 3, ... m 
^•- U 7=l,2, 3,...m+rc, 



except when i=j, in which case 



'h, i V 2j 2 • • • V m , m •*■• 



The unit-matrix \\v\\ therefore arises from the composition of two unit-matrices of the 
forms 

1 ,0 ,0 , . . , 0, 0, . . 

,1 ,0 , . . , 0, 0, . . 

,0 ,1 , . . , 0, 0, . . 

"•2, 19 "•2, 25 ^-2, 35 • • ^-2, ml "j ■*-, . . U 
"■», 15 "«, 25 "m, 3? • • ^B, m V, U, t • 1 



and 
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1, 0, 0, .. ,0 , ..0 
0, 1, 0, . . ,0 , . . 
0, 0, 1, . . ,0 , . . 
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0, 0, 0, .. u ltl , « 1>2 , ..u l>n 



0, 0, 0, . . u % „ w 2> , 



, u 



2, n 



"j «j U, . . u ni j , u„ i2 , . . u n>n 



||K| denoting a matrix of the type n X m of which the constituents may be any numbers 
whatever, and \u\ a unit-matrix of the type nXn. If for |]x|| we substitute the matrix 

HI=HI _1 xW, 

it is readily seen that we may invert the order of the factors in the expression of \\v\\ ; so 
that, using an abbreviated notation, the signification of which is evident, we may write 
either 



or 









1,0 

M 

1,0, 
0, u, 



X 



X 



1,0 

0, u 

1,0 

p, 1 



Substituting the latter expression of \v\ in the equation 



we immediately infer 



= Nx 



|4=||w|x||7+2foC||. 

c 

Every matrix satisfying the equation =1 is therefore comprised in the formula (55.) ; 

and since it is evident, conversely, that every matrix comprised in (55.) satisfies the 
equation, that formula contains the complete solution of the question. 

A particular solution of the problem (which may be taken for |y|) can be obtained 
as follows : — Complete the matrix |C|| by any n horizontal rows of constituents which do 
not cause the determinant of the resulting matrix to vanish. From this matrix a prime 
(i. e. a unit) matrix of the same type is to be deduced by the method of art. 3, a reduc- 
tion which can always be effected without changing the prime matrix |C||. 

Art. 11. The consideration of sets of fundamental solutions of linear systems is also 
of use in the theory of indeterminate systems containing terms not affected by any inde- 
terminate. Let 

A i) o+A ! . 1 # 1 +A t . 2 # 2 +...+A,. w+m =01 _ 

*=1, 2, 3, . . . n J 

represent such a system ; its general solution will assume the form 

jfc=l, 2, 3, . . . n+m) (57.) 

MDCCCLXI. 2 tr 



(56.) 
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where fflj, # 2 , ... a n+m is a particular solution of (56.), (&„ ^ 2 , ... p m indeterminate numbers, 
and ||a|| a set of fundamental solutions of the system 

-A-t, 1 ^T^i^ST • • • • T"A«, n+m^n+m == ^ \ /Kg \ 

t=l, 2, 3, ... to J ' 

Whenever, therefore, the proposed system is resoluble, its complete solution involves 
m indeterminates ; but in order that it should be resoluble, a certain condition must 
be satisfied by its coefficients. This condition is, " that the greatest divisors of its 
augmented and unaugmented matrices must be equal*." We shall call these divisors D 
and D respectively, representing the matrices themselves by |A|| and ||A |. That the 
condition is necessary may be seen by eliminating in turn every combination of to— 1 
indeterminates from (56). We thus find that every determinant of ||A|| is divisible by 
D , i. e. that D is divisible by D ; but evidently D divides D , so that D =D. To show 
that the condition is sufficient, as well as necessary, consider the system 



and let 



(59.) 
i=l, 2, 3, ... to 
(m+1) X(to+to+1) 



fl 



represent a set of its fundamental solutions. To say that 

(56.) is resoluble, is the same thing as to say that (59.) admits of solutions in which the 
value of w is unity; and (59.) will not, or will admit of such solutions according as 
9 lf , # 2i , .... m> do, or do not admit of any common divisor beside unity. But, by the 
theorem of art. 7, those determinants of ||d| into which the column 6 lt „, 2; . . . enters, 
are , equal to the determinants of ||A ||, taken in a proper order and divided by D. If 
D=D , the determinants of ||A„||, divided by D, are relatively prime, and consequently 
those determinants of ||0|| which contain $ h , Q 3y , . . . Q m< are also relatively prime ; a con- 
clusion which implies that d h0 , 8 2>0 , ... Q m0 are themselves relatively prime, i. e. that the 
system (56.) is resoluble. 

This criterion is not immediately applicable if the system (56.) be not independent, 
i. e. if the determinants of its augmented matrix |A|| be all equal to zero. But it may 

* [This Theorem has already been given by M. Ignaz Hegbe (Memoirs of the Vienna Academy, vol. xiv. 
second part, p. 111). I regret that in the abstract of the present paper, which has been inserted in the 
' Proceedings of the Eoyal Society,' no reference was made to M. Hegeb's Memoir, with the contents of 
which I was unacquainted, at the time at which that abstract was prepared. M. He gee's demonstration 
(adapted to the terminology here employed) is, in the main, as follows. (1) If the unaugmented matrix of 
an indeterminate system be prime, the system is always resoluble. For every determinate system, of which 
the matrix is a unit-matrix, is resoluble in integral numbers ; and we may suppose the given indeterminate 
system to form part of such a determinate system (see art. 10, supra). (2) The equation ||A||=||D|| x||A'||, 
in which ||D|| is a square matrix, having D for its determinant, and ||A'|| a prime matrix of the same type as 
||A||, is always resoluble (see art. 3). We can therefore replace the given system (56.) by a system of 
which the augmented matrix is ||A'||, and which is resoluble or irresoluble at the same time with the given 
system. But if D =D, the unaugmented matrix of this derived system is prime; i. e. if D =D, the pro- 
posed system is resoluble. (3) That the condition is necessary as well as sufficient may be proved as in the 
text.— Sept. 1861, H. J. S. 8.] 
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be applied to any independent system, equivalent to the proposed system, and deduced 
linearly from it. 

If we represent by Dj, the greatest divisor of the matrix, deduced from the matrix of 
(59.) by omitting from it the column A, i4 , A 2i , . . . A n>k , we may enunciate the follow- 
ing proposition : — 

" In every solution of the system (59.), the value of x h is divisible by -jy; and, con- 
versely, a solution of that system can always be assigned in which x h shall have any 
given value divisible by -jr." 

It will be seen that the solution of (56.) depends, first, on the solution of (59.), and, 
secondly, on that of the indeterminate equation 

0!, 0^ + 02,0 #2+- •.+fl»+l,0*»+l = l ( 60 -) 

If we represent the values of the indeterminates in this equation as the determinants of 

the matrix 

i=l, 2, 3, ... m 

,7=1, 2, ... ot+1 

(see art. 10), we may express the most general values of the indeterminates which 
satisfy (56.) in the determinantal form 

«1, ft «2, k "m + \,h 

«':~ 7i,i+f*i 9j,o> 7i,2+f*l 02,0, 7l,m+l +^l 0W+1.O 

72,1 + f*2 01,0? 72,2+f*2 02,05 7l,m + l + ^2 »m+l,« (^l-) 



7m,l+^)»0i,O5 7m, 2+f*m 02,0? 7m, ra+1 +/*»» 0«+ 1, 

Art. 12. We shall now indicate an important transformation of which any square 
matrix of integral numbers is susceptible. We begin with the following theorem : — 

" If a given rectangular matrix be premultiplied by a unit matrix, the greatest com- 
mon divisor of any vertical column of minor determinants is the same in the resulting 
as in the given matrix." 

For it is evident that any minor, either in the given or in the resulting matrix, is an 
integral and linear function of the minors formed from the same vertical columns in the 
other matrix. 

Similarly, it may be shown that 

"When a square matrix is post-multiplied by any prime rectangular matrix, the 
greatest common divisor of any horizontal row of minors is the same in the resulting 
rectangular matrix as in the given square matrix." 

For if 



nx(n-\-m) 

A 



nXn 
B 

2 u2 



X 



«X(i+ffl) 

c 



3L2 
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where ||C|| is a prime matrix, it is clear that every minor of |AJ| is a linear function of the 
minors formed from the same horizontal rows of ||B|| ; so that if a and b be the greatest 
common divisors of any corresponding horizontal rows of minors in those two matrices, 
a is divisible by b. But again, if Q be any one of the determinants of ||C||, and s be the 
order of the minors under consideration, any minor of ||B|, after multiplication by s , 
may be expressed as a linear function of a certain group of the minors taken from the 
same horizontal rows of ||A|. Consequently s X b is divisible by a ; or, since $ may have 
any one of a series of values which are relatively prime, b is divisible by a, i. e. b=a. 

By combining these results we obtain the theorem. 

" If v» Vm-h V«-2> • •• Vi represent the greatest common divisors of all the minors of 
order n, n— 1, . . . 1, respectively which can be formed out of a given square matrix, these 
numbers will remain unchanged, when the given matrix is premultiplied by any unit- 
matrix, and post-multiplied by any prime matrix whatsoever." 

!«x»! 



Art. 13. Let Q, the determinant of the square matrix 



, be a positive number, 



different from zero. It may be shown that by post-multiplication with a properly 
assumed unit ||«||, the matrix [a] can be reduced to the form 



(62.) 



/*I5 ^1,25 T\, 3 • • 


•• r,,„ 


5 fa 1 r 2, 3 • • 


• • ^2, n 


0,0 , p s 


• • ^3, n 


0,0,0 . . 


•• f>n 



where ^ 15 y> 2 , ... \h n are positive numbers, such that ^X^X-- X[6 n =Q, and the consti- 
tuents r iy k satisfy the inequalities 

0<r, ft <^ (63.) 

This was first observed by Gauss for the case n=2; by Sbeber for n=S; and the 
general theorem has been enunciated by M. Hbbmite*. Its precise statement is 

" Every matrix of the type nxn is equivalent (by post-multiplication) to one, and only 
one, of the reduced matrices included in the formula (62.)." 

To show this, let v h „ v % i,-..v n>l be the integral and relatively prime numbers which 
satisfy the equations 

fc2,3,...» )' <"■) 

and the inequality 

«!,» fll, 1+01,25 «2,1+ ■••«!,», ««,1>0. 



Then it is evident that, if \\v\\ be a unit-matrix of which «,,„ v 2il , ...v n>l form the first 
column, the matrix lldl X IMI will assume the form 



* Gauss, Disq. Arith. art. 213; Sbebee, " Untersuchungen ueber die Eigenschaften der positiven 
terniiren quadratischen Eormen" (Mannheim, 1831), art. 31 ; M. Heemite, Ceelle, vol. xli. p. 192. 
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^15 #1, 25 "I, 35 • • • V\, n 
« 5 02, 25 "2, 35 • • • 02, M 



^ 5 03, 25 03, 35 • • • 03, n. 



^ 5 0H, 25 "re, 35 •••0m 



where f*,=a lt „ «,, ,,+a,,,, %,, , + ...«,,„ ^,i- 
If this matrix be post-multiplied by the unit, 

1 , #25 "^35 • • • H>n 

0, 1, 0,...0 
0, 0, 1,...0 
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(65.) 



0, 0, 0....1 



the constituents b hi will be changed into h^^^Jc,, while all the other constituents will 
remain unaltered; so that by assigning proper. values to the numbers k 2 ...k n , we may 
bring the given matrix ||«|| into the form 

f*u T lt 2 , T lt 3 , . . . V 1> n 

*-> 5 02, 25 02, 35 • • • 02, re 
" 5 03, 25 03, 35 • • ■ 03, re 



" 5 0», 25 0», 35 • • • 0«, re 



where r^ ( verifies the inequality 



0<r, 



'pi 



From this it is easy to infer that if a matrix of the type (n— l)x(w— 1) can be reduced 
to the form (62.), the same reduction is possible for a matrix of the type nxn, i. e. since 
that reduction is possible when n=l, ?^=2, .. it is possible for every value of n. 

To prove that |]#| is equivalent (by post-multiplication) to only one of the reduced 
matrices (62.), it is sufficient to show that no two reduced matrices can be equivalent. If 
||a|| and |«'|| be two reduced matrices, and \\v\\ a unit-matrix, such that ||#||x||fl||=||ffl'|, it may 
be inferred, by comparing the corresponding constituents of the two matrices ||«|| X \\v\\ and 
\\a'\\ (beginning with the lowest horizontal rows of each and proceeding upwards), that 
all the constituents of |t>|| which lie below its principal diameter are zero ; and conse- 
quently that the constituents of the principal diameter itself are all positive units. 
Further, that the constituents above the principal diameter of ||«|| are likewise zero, may 
be established (for each line of constituents parallel to the diameter, beginning with 
that nearest to it) by means of the inequalities (63.) which are satisfied by the consti- 
tuents both of \\a\\ and \\a'\\. It thus appears that two reduced matrices cannot be equi- 
valent, without being identical. It will be observed that the reducing unit is unique ; 
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i. e. that only one post-multiplying unit can be assigned by which a given matrix can be 
reduced to the form (62.). 

If instead of reducing the given matrix \a\ by post-multiplication we employ a pre- 
multiplying unit, we obtain the following theorem : — 

" Every matrix of the type n X n and of determinant 6 is equivalent (by pre-multipli- 
cation) to one, and only one of the matrices included in the formula (62.), in which 
(«.„ \Jj % . . are positive, /«.„ (i 2 . . ia„=$, and r i>k satisfies the inequality 

0<r, 4 <^ 4 ." (66.) 

Art. 14*. The transformation to which we have referred in art. 12 is obtained by 
employing simultaneously a pre-multiplying and a post-multiplying unit-matrix. It is 
expressed by the equation 

|U||_ |U|v UZs-, V"-' , Vi vllfill (fi7 1 

||C*| — ||W|| A I » > . .. — A]|P||, V ' V 



_, V"-i , Vi 

V»-I Vra-2 Vo 

in which ||«|| is a given square matrix of the type »x», |«|| and |||3| are unit-matrices, and 
V«» V»-n Vn-si • • • Vd Vo are the determinant and greatest common divisors of the 
minor determinants of |«||, so that, in particular, v« is the determinant of |«||, y„_, the 
greatest common divisor of its minor determinants of order n— 1, Vi the greatest 
common divisor of its constituents, and Vi^l- The units ||aj| and ||/3|| are not abso- 
lutely determined, but admit, when n>l, of an infinite number of different values. 
H n=l, it is evident that the formula (67.) is verified; for we have the identical equa- 
tion (|a||=||l||x — X(|l||. It is therefore sufficient to show that, if the transformation 

indicated in the formula can be effected for matrices of the type (n— l)x(w— 1), it 
can also be effected for matrices of the type »x». The demonstration depends on an 
elementary principle, which it is worth while to enunciate separately. 

"If 

Ui=A,- (1 ^i+A,- i2 a? 2 + . . . -f- A t - ; n + m %„+ m , 1 

t=l, 2, 3, ...» J ( -) 

denote a system of n linear functions of »+m indeterminates, (m>0)„and if the consti- 
tuents of the matrix |A|| do not admit of any common divisor, it is always possible to 
assign integral values to oc x , x a , ... cc n+m , which shall render U„ U s , . . . U„ relatively 
prime." 

For, in the first place, we can obtain values for U„ U 2 , . . . U„, which shall not have 
any common divisor with a given number M. Let p, q, r . . . be the different prime 
divisors of M; one at least of the constituents of ||A||, for example A iiP is prime top. 
Attributing to Xj a value prime top, and values divisible by^> to the remaining indeter- 
minates, we shall obtain for U 4 a value which is certainly prime to p. Similarly, by sub- 
jecting the indeterminates to proper congruential conditions with respect to the modules 
q, r, . . ., we can render one, at least, of the functions U prime to q, one prime to r, and 

* [This article has been in great part rewritten since the paper was read. The demonstration is not 
essentially changed, hut is presented in what seems to be a simpler form.— Sept. 1861, H. J. S. S.] 
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so on ; i. e. since we can assign to the indeterminates values simultaneously satisfying 
all these congruential conditions, we can give to Uj, U 2 , . . . U„ values the greatest 
common divisor of which is prime to M. Let D n be the greatest divisor of j|A|j, D H _, 
the greatest common divisor of the first minors of ||A|| ; and let C„ C 2 , . . . C n be a set of 
simultaneous values of U„ TJ 2 , . . . U„, having a greatest common divisor c, which is prime 

to ty"— . Since the equations 

^,AT-^f,AT •• • "T"A' ( n+mfin+m-— ^is 

i—1, 2, 3, ... n 

are resoluble, it will follow from the condition of resolubility (see art. 11), that the 
determinants of its augmented matrix, and in particular those which contain the column 
C„ C 2 , . . . C„, are divisible by D n . Let ^XcXD„_, be the greatest common divisor of 

these last determinants; then ^x<?xD„_, is divisible by D M , i. e. 6 is divisible by t~. 
It appears from this, that the condition of resolubility is satisfied by the system 

C- 

A', 1#I "T-"-f, 2^2 "T" •• • T-"i ) »+m'''ii+m = 7' 

i=l, 2, 3, ... n, 
that is to say, it is possible to obtain a simultaneous system of relatively prime values 
forU„U 2 , ...XL, 

To apply this principle to the transformation of the matrix ||<zj|, let 

fo'l^g 2 . W 

Vn-I aa i,j 

The constituents of the matrix || [a] || do not admit of any common divisor ; consequently, 
in the system 

K ,]*,, .+[« 2 , 0^2, i+ • • • +[a n , $», i=% „ \ 

i=l, 2, 3, ...» /' 

we can assign values to b h „ b % „ . . . b n< n which shall render u h x , u 2) ,, . . . u n> j relatively 
prime. Let ||w|| denote a unit-matrix of which the first column is w,_ „ w 2> „ . . . u, h , ; 
and p|) a square matrix of which the first column is b h „ 5 2i „ . . . b n> „ and of which the 
remaining constituents are defined by the equations 

»=1, 2, 3, ...» y (71.) 

j= 2, 3, ... n. J 

Observing that the systems (69.) and (70.) involve the inverse system, 

a i, I w l v l T#», 2^2, j T • • • T fl t, m**m, 1 == y "t, » I f72 1 

e=l, 2, 3, . . . n, J 

we infer that the matrices j|w|| and \\b\\ verify the equation 



Hx|MI=#|x 



v * i i 

, X, i, . . . 



Vm-1 



(73.) 
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in which 
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V« 



V«-i 



,1,1,-.. 



denotes a matrix of the type nxn. It follows from (73.) 

that v»-i is the determinant of \\b\\ ; let that matrix be reduced by premultiplication 
with a unit-matrix ; and let 

IHI=Hx||v B -,||, (74.) 

where ||«|| is the reducing unit, and ||v»-i|| the reduced matrix, 

ft, **1,2, ^1,3 •• • T\,n 



" , ft , *"s, 3 • • • *"a, n 

0,0 , ft ...r 3>n 



0,0,0 
so that (73.) assumes the form 

H=IMIx||v-»|x 



. (75.) 



v» 



v»- 



ft 



,1,1,... 



xM 



(76.) 



It may be proved that in (75.) ft=l, r 1)2 =0, r it3 =Q ..r hn =0. For since the matrix 



V.-X 



— , -L, X.. 



is derived from ||«| by multiplication with unit-matrices, v«-i is the 



greatest common divisor of the first minors of |v»-i|| X 



. V» 1 1 



Therefore y„. 



divides ft X ft X • • X ft, which is one of those minors ; but also Vn-^ft X ft • • X ft ; «• «• 



V« 



V»-i 



,1,1,.. 



assumes the form 



^,=1, ftXftX..Xft=V ra _,, and the product ||v»-i||X 

n , '1,2, ' l.g, • • • ' l,n 
Vm-1 

V ) ft I'^U"' ^*2, » 

, , ft , ... r 3>n 

, , , . . . ft 

One of the minors of this matrix is r 1)2 Xft . • • Xft, which cannot be divisible by v»-i 
or ft X ft X ... X ft, unless r lt 2 is a multiple of ft; but r 1( 2 < ft, because ||v»-i|| is reduced, 
therefore r, )2 =0. Similarly, it may successively be shown that r h3 =0 . . . r 1>n =0. Now 
if the matrix 

ft, 7*2, 3, ' • • *2, 



, ft , . . . r, 



3, « 



(77.) 



0,0 , ... ft 
which is of the type (w— l)x(%— 1), be reduced to the form 



\w\\ x 



V«-I V«-2 
V»>-2 Vn-3 



Vo 



X 



w 
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in which v»-2, V«-3, • • • represent the greatest common divisors of the minors of (77.), 
we may replace ||v»-i| by the matrix 



V»-.= 



10 
0t/ 



X 



1, 



V»-i V»-a 



V 

Vo 



X 



10 
0«' 



where 



10 

Ov 1 



and 



10| 

0u' 



Vra— 2 V«-3 

denote unit -matrices of the type nxn, the forms of which are 



sufficiently indicated by the symbols themselves. Hence, observing that 



and that 



10 
0w' 



1, 



XI 

I 

V»-l Vm-2 



Vm- 



, 1, 1, 



V»-i 



,1,1,... 



X 



10 
0u' 



V 
Vo 



V«-2 V»— 8 

we obtain, from (76.), 

H=.IMIx 

or more simply, 

IHI=Nlx 



x 



v« 



V»-i 



,1,1, 



Vm Vre-1 



10 

0v' 



X 



V» Vw-i 
V«-i' V»-s' 

V» V«-l 



V»-l V«-2 



V 
Vo 

V 
Vo 



X 



V»-l Vm-2 

0«' Xl[u " 



V 
Vo 



x 



It has, however, still to be shown that v»- 2 , V»-s, • • • which have been defined with 
reference to the matrix (77.) are the greatest common divisors of the successive systems 
of minors of ||a||. These greatest common divisors are the same for the given matrix |«| 

which is derived from it by multiplication with 



and for the matrix 



Vn Vn- 



Vo 1 



V»-i V»-2 

unit-matrices; consequently v»-i divides every first minor of 



V» Vn-l 



V»-l Vn-2 



V 

Vo 



, and, 



in particular, it divides -^ X^X^X... X^= V -~^=* ', i. e. ^ divides -^- 

V»-l Vn-S Vn-4 Vo Vn-l V»-2 V»-i' 

Again, V»-n Vn-a> • • • Vi, Vo, which are the determinant and greatest common divisors 
of the minors of (77.), are also the determinant and greatest common divisors of the 
minors of the matrix 

y»-i v»-2 Vi 

Vn-2 Vn-3 Vo 

so that if s< n — 2, y„ divides every minor of order s in (78.), and, consequently, the minor 
— — X ^^ X zr- 2 X .... X ^r ; or ~- divides -—-. It thus appears that in the series of 

V* V«-2 Vs-3 Vs V»-i V« 

numbers 



(78.) 



V«-i 



V? Vi 
Vi' Vo 



V»-i V»-s 

each term is divisible by that which comes after it. Every product of s terms of that 
series is therefore divisible by the product 



V, 



zr J -X-..X^r=Vs'> or, which is the 



V»-i Vs-2 " Vo 

same thing, v s is the greatest common divisor of the minors of order s in the reduced 
matrix (78.), and therefore in the given matrix ||«||. 
mdccclxi. 2 x 



518 
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Art. 15. If the proposed matrix |«|| be not square, but of the type nX (n-\-m), let 
IHI = ||v*|XJ|«'|| 5 where ||a'|| is a prime matrix of the same type as ||a|, and |v„| a square 
matrix, of which the determinant is v»> the greatest divisor of ||<z||. Then if ||v»|| be 
expressed in the form 



NX 



Vn Vn-1 



Vj 
Vo 



xW 



I Vra— 1 V«— 2 

and if, for brevity, we write |V| for ||w| x |a'|[, we obtain for |«|| the expression 



W=iX 



V» V«-i 



Yj 

Vo 



X||V| 



r v „ .. .. (79.) 

V«-l Vn— 2 

The numbers v»> V»-» • • ■■> which are the greatest common divisors of the minors of 
|vj, are also by the theorem of art. 12, the greatest common divisors of the minors 

of \\a\\. We see therefore that -Yi. is always divisible by Yt=i t in the case of an oblong 

y*-i Vs-2 

as well as a square matrix. 

Art. 16. To show still more clearly the nature of the quotients -5i-, we add the 

Vs-1 

following proposition : — 

" If in any rectangular matrix we divide each minor determinant of order s by the 
greatest common divisor of its own first minors, the greatest common divisor of all the 

quotients thus obtained is ^r J --" 

By this proposition — — is itself defined as a greatest common divisor, instead of being 
defined as the quotient of one greatest common divisor, divided by another. 



To establish its truth Ave may first consider the quotient 



Vn 

V»-i 



in any rectangular 



matrix j|Aj| of the type nx(m-\-n). Let u denote the greatest common divisor of the 
quotients obtained by dividing each determinant of |A| by the greatest common divisor 
of the first minors of that determinant : we have then to show that 



Vn 



Vn- 



■ = &). 



Since the greatest common divisor of any vertical column of minors in ||A|| is not 
altered by premultiplication with a unit-matrix, it is evident that w as well as -^- 2 - 

Vn-I 

will remain unchanged by that operation. ' If, therefore, 



NX 



Vre Vra-1 
, , 



y 

Vo 



xllVII 



(79.) 



V»-l Vra-2 

where \\ is a unit, and ||V|| a prime matrix, we may consider instead of ||Aj|, the simpler 
matrix 

Vn Vra-1 Vl' 



X]|V||. 



(80.) 



Let ll^ll, 



Vn-l Vn-2 Vo 

&c. be the different square matrices of ||V|; n 2 , . . . their determi- 
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nants ; tpi t ne greatest common divisor of those first minors in |0J which do not contain 

the constituents of its uppermost row, so that if is integral ; lastly, let w i be the quotient 
obtained by dividing the determinant 



Vra Vm-1 



Vl 

Vo 



xpl 



v _ „ (81.) 

Vm-1 Vm~2 

by the greatest common divisor of its first minors, so that a is the greatest common 

divisor of a lt u 2 Now the greatest common divisor of the first minors of (81.) is 

evidently divisible by v»-» an< i divides y»-i X ^ (because Vn-i^i is the greatest common 
divisor of one of its rows of minors). Consequently w ; divides V»A-J-V»-» an d is divisible 

hy vJi + Vn-i^i- Therefore — — is a common divisor of certain numbers respectively 



Vm-1 



Vm &l 



dividing the numbers «„ » 2 , . . ., viz. the numbers — — •-/-; it is also (because 0„ 6 2 



are 



V» 



relatively prime) the greatest common divisor of the numbers — - { , in which the same 



Vm-1 



Vn 



numbers <y t - are respectively contained ; i. e. — — is the greatest common divisor of the 

Vm— i 

numbers <y i5 <y 2 . . . themselves, or 



-■=■&>. 



Vm-1 

By the aid of this particular case of the theorem the general proposition itself may be 
proved as follows : — 

If in any rectangular matrix of the type nX(m-\-n) we propose to determine Cl„ the 
greatest common divisor of the quotients obtained by dividing each minor determinant 
of order s, by the greatest common divisor of its own first minors, we may begin by 
selecting any s vertical columns [s<n], and forming the proper quotient for each deter- 
minant of order s, contained in this partial matrix of the type*w x s. Let K t denote the 
greatest common divisor of these quotients ; then, as we have just seen, Xj is the greatest 
common divisor of all the determinants of the partial matrix, divided by the greatest 
common divisor of all its first minors. Hence (by art. 12) X; will remain unchanged 
when the given matrix is premultiplied by a unit-matrix. But Cl s is the greatest 
common divisor of all the divisors X 15 X 2 . . . corresponding to every group of s vertical 
columns ; therefore Cl $ is itself unchanged by premultiplication. Similarly, if a square 
matrix be post-multiplied by a rectangular prime matrix, it may be shown that Cl s is 
the same for the given square matrix, and for the resulting rectangular matrix. Hence 
if, as before, 

Vn Vl 

Vtt-i Vo 



IIA||=Hx 



£1, and -¥±. are the same for 



Vm 

|V»— i 



Vs-1 

Vo 



V« V«.-i 

|Vn-l Vm-2 



V 
Vo 



X||V||, 

and for |A|. 



But in the matrix 



it is evident that -¥±- and Cl s coincide ; therefore in any rectangular matrix 

Vs-J 

V.-i s 

2x2 
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From the definition of — — as a greatest common divisor, which we have now obtained, 

we infer that if ||D| be any matrix containing another matrix ||y|, and if D s , D s _! . . . v« 
V,_„ ... be the greatest common divisors of the corresponding minors in ||D| and |v|| 

respectively, not only is y s divisible by D„ and v,_, by D,_„ but also -^- by =— — • 

It is not difficult to show that in any matrix -^ £ - is the greatest common divisor of 

all the quotients obtained by dividing each minor of order s by the greatest common 
divisor of its minors of order s— Jc. But as this extension of the preceding result is not 
needed in what follows, we may omit it here. 

We may add, that the theorem of this article is precisely equivalent to the following, 
which may be demonstrated by a different method. 

"If P * be the highest power of a given prime that divides all the minors of order s 
in a given matrix, and if all the minors of order s— 1 contained in one particular minor 
of order s are divisible by P s_1 , that minor is itself divisible by P s ." 

It should be observed that whenever all the minors of any determinant are zero, the 
quotient obtained by dividing the determinant by the greatest common divisor of its 
minors is also zero. 

Art. 17. These results admit of immediate application to the theory of systems of 
linear congruences. The general type of such systems is 

A itl w l +A ii ,sc a +.. .A i>n x n =A i>n+ „ mod. Ml 

i=l,2,3,...n> J ; (82 ° 

and to construct a complete theory of them it is requisite, first, to assign a criterion for 
their resolubility or irresolubility ; secondly, when they are resoluble, to investigate the 
number of incongruous solutions of which they are susceptible ; and, lastly, to exhibit a 
method for obtaining all these solutions. "We shall first suppose that n'=n ; i. e. that 
the proposed system is neither defective nor redundant. 

Let D ra , D„_j . . . . v« Vn-u • • • respectively denote the greatest common divisors of 
the determinants and minors of the augmented and unaugmented matrices of the 

system (82.); also let i„, &„_„ . . . &, denote the greatest common divisors of M with -^-, 

Vn-l 

of Mwith zr- 1 , ■■■, and let d n , d n _ x ... similarly represent the greatest common divisors of M 

\n— 2 

with f~, of M with ff=*., Sec. ; then, if d=d„X d n _, X . . . . X i„ ^=Kx\-, X . • . X K 

we have the two following theorems : 

(i.) " The necessary and sufficient condition for the resolubility of the system (81.) is 

d=r' 

(ii.) " When this condition is satisfied, the number of its incongruous solutions is d." 
To demonstrate the first of these theorems, we revert to the principle of art. 11, 

from which it appears that the necessary and sufficient condition for the resolubility of 

the system (82.) is that the greatest divisors of the two matrices 
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and 



M,0 J 0,.0,A liI ,...A l> . 

0,M, 0,.0, A 3>1 , ...A 3> „ 
0,0, M, . 0,A 3>1 , .... A,, fl 



0,0,0, 


. M, A.,,, . 


• • A„ ; „ 


M, , , 


0,A liI( . 


• • A, n+l 


, M, , 


, A 2) l5 . 


• • -^S.n + l 


, , M, 


» A 3> „ . 


• • Al, »+l 



(83.) 



(84.) 



0,0,0, .M, A,,,, ... A, >re+I 

are to be equal to one another. Now the first of those greatest common divisors is 
evidently the greatest common divisor of 

M", M- 1 v» M"- 2 v» • • • • M v»-n v. ; 

which, for brevity, we shall represent by the symbol 

[M",M- I Vi,M- 1 v»..-V.-,M,vJ (85.) 

Let M=PxQxE • • ., P, Q, R, • • • denoting powers of different primes ; we may then, 
in (85.), replace M by P, Q, P, . . . successively, since 

[M-.M—Vm-.-Mv-lV.] 

= [P», P"- 1 V „ P V.-» V.] X [_Q n , Q"" 1 V, • • • V„] X • • • 

If P divide any one of the numbers — — . . . . , let -~- be the least of them that it divides ; 
J V«-i V»-i 



also let Pj= 



P, 



V«-i 



; so that P(=P, if i>s. Then 

[P» ! F- 1 V„....Pv»- 1 ,Vn] 

P s/1— pn-lVl p«-2V2 V»l 

— riX |jv Pi' p 1 '--p 1 J 



= P:X 



p»-l p«-2 Vg pn-3 V3 

r ' r V,' X Vi' 






observing that ty- 1 is prime to P [if s>l], and that we may therefore divide the last n 

_ p 

numbers by p- 1 ; and may then omit p^, which is divisible by P M_1 . Continuing this 

process, we find 

[P",P- , v.,.-..Pv.-»v.] 
=p, xP 2 x . • . xP.- 1 [p— + \ P"-^, P — 1 g £ii , • • • ^- 

(_ Vs-i Vs-i V»-i 
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or, since £- is divisible by P, and 2j±* =-2tt*. x Vf±*n . . . x Jt- bj p* + ., 

V.-i ' V«-i V*+*-i V*+*-a V«-i 

CP»,P"- , v.,.P—v»....Pv^i,v.] 
=p 1 xP 2 xP 3 ... .xP^,xp— +1 
— n'~" p 

But ^=P i xQ,XRjX • ••; and consequently the greatest common divisor of the deter- 
minants of (83.) is fcjX^x . • • X&„ or &• Similarly, the greatest divisor of (84.) is 
diXd 2 X . • . Xd„ or d. The necessary and sufficient condition for the insolubility of 
the proposed system of congruences is therefore contained in the formula 

d=l. 

It should, however, be observed that, since =p— - divides -^~ (art. 16), d s divides 5 S , and 

therefore the equation 

d=l 

involves the coexistence of the n equations 

d x —\, d 2 =h 2 , . . . d n =\. ........ (86.) 

To investigate the number of solutions of the system (82.), supposed to be resoluble, 
let ||«|| and ||/3|| be two unit-matrices satisfying the equation 



Vra V»-l V- 



V«-i V«-2 Vo 



(87.) 



Hx||A||xP||= 

also let 

**= A, i »i +ft, ■«.+■• ■ ft, . «. ) 
»=1, 2, 3,...» j 

j=l, 2, 3, ...%. j 

Then it is evident that the proposed system of congruences is precisely equivalent to the 
system 



Vn~i+ 



-Vi—iCi, mod. M, 



t=l, 2, 3, ...» J 

in such a manner that the two systems are simultaneously resoluble or irresoluble ; and 
that from any number of incongruous solutions of the one an equal number of incon- 
gruous solutions of the other is deducible. But the whole number of incongruous solu- 
tions of (88.) is S, x<i a X-..X&„=&; i- #• the number of solutions of the proposed system 
is S. 

By the use of the unit-matrices ||aj| and |[/3|| the actual resolution of the proposed 
system is made to depend on the resolution of the n congruences contained in (88.). 
But this method of solving a system of linear congruences, though very symmetrical, 
is perhaps too tedious for the purposes of computation. 
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Art. 18*. Let the proposed system of congruences be the defective system 
A,. ;I ar 1 +A ei2 ^ 2 +...+A I)B+m ^ +m =A ii „ +m+1 , mod; M,l 

i=l, 2, 8,...», J ( '' 

and let the notation of the last Article be retained. It is easily seen that the condition 
of resolubility of the system (89.) is, as before, 

l=d. 
But the number of its incongruous solutions, when that condition is satisfied, is not l, 
but 5 X M m . For we have seen that we can find a unit-matrix ||«||, and a prime matrix 
|A'|| of the type nx(n+m), satisfying the equation 



WX||A||= 



v« 


v»-> 


Va Vi 


V»-i' 


V»-a' ' 


' ' Vi' Vo 



X || A||; 



^ , U < =C,,mod.M, (90.) 

Vn—i 



we may therefore replace the system (89.) by a system of the form 

in which 
and 

If the system (89.) is resoluble, the system (90.) will be so too, and will give d or & different 

systems of values for TJ„ U 2 , . . . U„, any one of which may be represented by the formula 

V ( =Ui, mod. M, 1 

I (91.) 

Let us replace the modulus M by P, the highest power of one of its prime divisors. 
Since IJA'I is a prime matrix, one at least of its determinants, for example, the deter- 
minant 2+ A',, j A' 2) 2 . . . A' nj „, is prime to P. It will follow from this that, whatever values 
we attribute to x n+l , x n+% , . . . w n+m , each of the S systems represented by (91.) is resoluble 
for the modulus P, and gives, for any assumed values of x„ +l , x n+2 , . . . x n+m , only one set 
of values of #„ x 2 , . .. x n . Each of those I systems admits, therefore, of P m solutions for the 
modulus P, i. e. of M m for the modulus M. The system (89.) will consequently admit 
of 5xM m solutions. 

Let us also consider the redundant system of congruences, 

A ! . jl # 1 +A i(2 #2+....+A i)n # s H=A,. mod. M,] 

y .... (92.) 
•=1,2, 3,...w+m, J v ' 

and let D„ + , denote the greatest divisor of its augmented matrix. Let p represent a 
prime divisor of M, and let p e , p\ p is be the highest powers of p, which divide M, D,, V« 
respectively. The condition of resolubility of art. 11, applied to the system (92.), con- 
sidered with respect to the modulus p e , becomes, after division by p (m ~ l)e , 

[^ I »+i,^» I » +8 ,j? I «-i +2e , j?(»+ J )9]1 

=[j) i » +9 ,/»-i +M , p(»+i)ejj' •••••• i - •; 

* [This article has heen added since the paper was read. The theorems contained in it are supplementary 
to that of the preceding article. September 1861, H. J. S. S.] 
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And this equation is impossible, if 0>I„ +1 — 1„. For I s+l — I S >I S — I,_„ because -f^ 1 is 

divisible by fp 1 -; the inequality, I n+ i<I«+^ involves, therefore, the inequalities 

I n+l <I n -s+i+S0, 1 

«=l,2,8,...»+lj 

and these, again, imply the corresponding inequalities 



(95.) 
s=l,2, 3,...»+l,J 

because I„_ s+1 < V-s+i- From (94.) it appears that the value of [j> 1 » +1 , p 1 "* 9 , • • .p (n+1)e '] 
is <p l »+i, and from (95.), that the value of [p l " +e , jA»-> +29 , . . ._p (re+1)< '] is a power of p 
superior to^A+i ; i. e. the equation (93.) is impossible. We thus obtain, as a first con- 
dition for the resolubility of the proposed system (92.), the congruence 

^^O, mod. M (96.) 

When this condition is satisfied, we obtain from (93.), omitting the term jpWi (because 
I„+^I»+^)j an d dividing by j) 9 , the equation of condition, 

= (j»S y»-i +9 , y»-2+ 29 , . . . y 9 ], 
which leads us (as in the last article) to the simple formula 

This equation, therefore, and the congruence (96.), express the necessary and sufficient 
conditions for the resolubility of the proposed redundant system. 

When these two conditions are simultaneously satisfied, the number of incongruous 
solutions is I. For, if we again consider the proposed system of congruences with 
respect to the modulus p 9 , and select from it a partial system of n congruences such that 
the determinants of its augmented matrix, which are necessarily divisible by p\ are not 
divisible by any higher power of p, it is readily seen that every set of values of the 
indeterminates #„ x 2 , . . . x n , which satisfies the partial system, will also (by virtue of 
the inequality ^<I M+1 — 1„) satisfy the remaining congruences of the proposed system. 
The number of solutions of the proposed system is therefore the same as that of the 
partial system. And because p 1 " the highest power of p which divides every determinant 
of order n in the augmented matrix of the proposed system is also the highest power of 
p which divides the augmented matrix of the partial system, it follows from the last 

theorem of art. 16, that p l «-i, p l "-% are the highest powers of p which divide the 

corresponding orders of determinants in the latter, as well as in the former matrix. The 
number of solutions of the partial system (and consequently of the proposed system), 
considered with respect to the modulus p 9 , is therefore expressed by the formula 

\_p ln , p l «-i+o, jp" 8 ]; 
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or, finally, the number of solutions of the proposed system, considered with respect to 
M as modulus, is d or 5. 

Art. 19. We shall terminate this paper with an elementary theorem relating to linear 
systems of equations, which admits of frequent application in other parts of the theory 
of numbers. 

Resuming the notation of art. 11, we may see from the theorem of that article, that 
if the system (56.) be resoluble for any given values of the numbers A 1)0 , A 20 , ... A Wj0 , it 
is also resoluble for any other values of those numbers, respectively congruous, for the 
modulus D, to the given values ; so that the resolubility or irresolubility of the system 
depends exclusively on the residues of the numbers A i; , mod. D. There are D M possible 
combinations of these residues, and we shall now show that for D n_l of them the system 
is resoluble, while for the remaining D" _1 (D— 1) it is irresoluble. For this purpose let 



|«||X|A||= 



Dn-l JL'm-2 *J, 



xfA'l, (97.) 



||«| denoting a unit-matrix, and ||A'|| a prime matrix of the same type as ||A||, while D M , 
D«-i5 Di, Do ar 6 of course the greatest common divisors of the determinants and minors 
of ||A|. Let also 



The given system is then exactly equivalent to the system 

i=n. 1 

. . . (98.) 



Dl i. n %i~T A. 2 &2\ • • •~V A, n+m <^n+m\ — ^t 
«— i 



«=1,2, 3,...w. 



J' 



For the resolubility of this system it is requisite that Q should be divisible by -^ t+1 • 

and this condition is sufficient as well as necessary, because ||A'|| is a prime matrix. Now 
of the D or D n values, incongruous for the modulus D, which may be attributed to C e , 

-t| — are divisible by T ^~' +1 -; whence it is evident that of the D" systems of values which 



may be attributed to C„ C 2 , . .. C„, D M -r- 



, i. e. D" -1 render the system 



Dj Dg D3 T> n 
D D 1 D 2 "D n - K 

(98.) resoluble. Consequently the given system is also resoluble for D n ~\ and no more, 

of the systems of values that can be attributed (mod. D) to A, s , A a> B , ... A m> . 

Art. 20. The methods employed in the present paper are without exception such as to 

be immediately applicable to any species of complex numbers which admit of resolution 

into actual or ideal prime factors. And the greater part of the results at which we have 

arrived may be transferred, mutatis mutandis, to the theories of such numbers. For 

example, if in the equations (56.) we suppose the constituents of ||A|| to represent complex 

numbers, it will be found that the criterion for the resolubility or irresolubility of the 

system, which we have demonstrated in the case of ordinary integers, applies equally in 

the case of complex numbers ; and again, the condition of resolubility of a system of 

congruences of which the modulus as well as the coefficients are complex numbers, is 

MDCCCLXI. 2 Y 
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precisely the same as in the case of common whole numbers ; while the expression for 
the number of the solutions (when the condition of resolubility is satisfied) is simply 
the norm of m. 

But without entering into the developments which this extension of the subject of 
this paper would require, we shall confine ourselves to an application of the result of 
the last article to a demonstration of the fundamental principle in the arithmetical 
theory of complex numbers, that the number of incongruous residues for any complex 
modulus is represented by the norm of the modulus. 

Let a, be one of the roots a„ a 2 , ... a n of the equation T? n (x)=0, which is supposed to 
be of n dimensions, to be irreducible, and to have all its coefficients integral, and that 
of its first term unity; Let also <p M _,(«) be the complex modulus under consideration; 
its norm, which we shall symbolize by N, is defined by the equation 

N=N.<p M _ ] («)=n::r<p B _. 1 («o. 

Consider the N 2M _, residues (incongruous mod. N) which are included in the formula 

E 2M _ 2 ( a ), ........... (99.) 

where Ea„_ 2 denotes an integer function of 2n— 2 dimensions; it is evident that every 
complex number is congruous, for the modulus <p„_,(a), to one at least of these N 2ra_1 
residues. If E and E' be any two (the same or different) of the same residues, it is 
also plain that the congruence 

Rs^R', mod. <P„_i(oj) 
will, or will not, be satisfied, according as it is, or is not, possible to assign two functions 
of os, F„_,(#) and <p re _ 2 (#) having integer coefficients, and satisfying the equation 

F.(ar^ a (*)+F.. 1 (*)^_ 1 (ar)=E(a;)-E'(ar) (100.) 

This equation is equivalent to a system of 2n—l linear equations, in which the 
unknown quantities are the 2n—l coefficients of <p M _ 2 (#) and F w _,(#), and of which the 
determinant is the dialytic resultant of F w (x) and <p M _,(#), i. e. the norm of <p„_i(«) or N. 
If then we suppose E(a) to represent any given residue included in the formula (99.), it 
will appear from the theorem of the last article that the equation (100.) is resoluble for 
N 2 * -2 different values ofE'(;r), i.e. that every complex number is congruous, for the 
modulus <p w _,(a), to precisely N 2 "- 2 of the N 2 ™ -1 residues contained in the formula (99.), 
or that the number of residues, incongruous mod. <p B _,(«), is precisely N. 

It is, however, proper to observe, that a complete demonstration of this important 
theorem has already been given by Professor Sylvbstee (see a paper signed "Lanavi- 
censis," in the 'Quarterly Journal of Pure and Applied Mathematics,' vol. iv. p. 94 
and- 124). 



